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Supersymmetry in the gauge sector could be realized as N=l or N=2 Supersymmetry, but the 
current LHC searches assume an N=l realization. In this paper we show that squarks could be as 
light as few hundreds of GeV for N=2. We also describe an experimental procedure to count the 
number of supersymmetries, i.e. to distinguish between N=l and N=2 supersymmetry, based on 
counting bins with different jet multiplicities and number of leptons. 
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I. INTRODUCTION 

Supersymmetry, also known as SUSY, is a very attrac- 
tive proposal for physics beyond the Standard Model. 
SUSY provides a mechanism to stabilize the electroweak 
scale, hence leading to a natural Higgs scenario. SUSY 
also offers an explanation for Dark Matter, and may even 
lead to unification, i.e., the idea that the three forces of 
the Standard Model are actually unified into a single force 
at high energies. 

SUSY's successes rely on new particles near the elec- 
troweak scale, particles which should be accessible at 
the Large Hadron Collider (LHC). The LHC is actively 
searching for SUSY, so far with negative results. Cur- 
rent bounds on superpartner masses already undermine 
the theoretical appeal of SUSY. 

But experimental searches make assumptions about 
how SUSY is realized at low energies. In particular, they 
are limited to only one SUSY generator, the so-called 
Af = 1 SUSY scenario, also known as the Minimal Su- 
persymmetric Standard Model (MSSM). More generally, 
the number of SUSYs is something we need to deter- 
mine. As we explain later, a very attractive possibility 
is that SUSY comes in a different version, called Af = 2 
SUSY pQ. 

In this paper we show that the current searches are 
very sensitive on the assumption of Af = 1 SUSY: al- 
lowing for Af = 2 SUSY opens a new area of the pa- 
rameter space where light supersymmetric particles are 
still allowed. We then describe a set of observables which 
can be used to characterize the SUSY scenario. In other 
words, we propose a procedure to count the number of 
supersymmetries. This characterization is based on the 
current searches of ATLAS and CMS of multijets plus 
missing energy, with and without leptons. 

To approach the problem of counting SUSYs, we ex- 
ploit the fact that Af = 1 SUSY implies that the gluino, 
the fermionic SUSY partner of the gluon, is a Majorana 
fermion, whereas Af — 2 SUSY leads to a Dirac gluino. 

The MSSM is the simplest version of SUSY, and is 
therefore traditionally adopted in most studies. However, 
there are many reasons to turn towards Af — 2 SUSY: it 
allows natural scenarios for Dark Matter [21 [3] without 
the well-tempered tension it provides interesting so- 



lutions to the flavor problem in SUSY [5 ; it has a softer 
behavior than Af = 1 in terms of UV structure [Sj ; gauge 
unification can be achieved i2, (i ;7j, among many inter- 
esting aspects of N = 2 SUSY [3]. 



II. MODEL BUILDING 

One can ask the question of how many SUSYs lead to 
a reasonable low-energy theory, and the short answer is 
that matter multiplets should come as Af = 1 SUSY, but 
Af = 1 or Af = 2 SUSY in the gauge multiplets are both 
viable possibilities, see Ref. [B] for a discussion on the 
subject. 

One could try to extend the matter sector to Af = 2 
SUSY, as long as the gauge sector couples non-linearly 
to the matter sector [B], but this is a contrived possibil- 
ity. One could also try to extend the gauge sector to 
Af = 4 or higher, but then the presence of so many extra 
fields give rise to Landau poles for the gauge couplings 
at an unreasonably low scale. One could also engineer 
by hand a new fermionic color octet with a Dirac mass 
term with the gluino, and with that reproduce the same 
phenomenology as described in this paper. This ad-hoc 
modeling would not have the advantages of Af = 2 SUSY, 
such as the supersoft behavior, or the Dark Matter fea- 
tures. 

For those reasons, we are going to compare Af — 1 and 
Af = 2 SUSY in the gauge sector, and our comparison 
will rely on the Majorana or Dirac nature of the colored 
gauge fermion, the gluino. 



III. GENERAL ARGUMENTS 

Unless there is a large gap between the SUSY colored 
and electroweak states, the strong production of squarks 
and gluinos is the dominant production mechanism for 
SUSY particles, and dominates over the electroweak pro- 
cesses, which we neglect in the rest of the paper. 

Strong production at leading order would lead to three 
different states 

• Squark pair-production (qq) 
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where q refers generically to a squark or an anti- 
squark q* . The squark pair production process is 
dominated by the i-channcl gluino exchange dia- 
gram, but in the case of same flavor squark anti- 
squark production also the s-channel gluon medi- 
ated process is available. If the gluino is a Dirac 
particle, it cannot mediate a chirality flip. There- 
fore the i-channel exchange of the Dirac gluino for- 
bids final states with squarks of the same helicity, 
and the production cross section of the qq final state 
is reduced in the Dirac case as compared to the 
Majorana case: only the final states with squarks 
of opposite helicities are allowed, i.e qLq' R or q* L (f L 
but not e.g. <?l<Zl- 

• Gluino pair production (gg) 

The gluino pair production process proceeds 
through the t-channel exchange of a squark (from a 
qq initial state) or a gluino (from a gg initial state) , 
or through the s-channel process mediated by a 
gluon (from a qq or gg initial state). Also here the 
t-channel diagram is dominant. The cross section 
of the gg final state is enhanced in the Dirac model 
as compared to the Majorana case, because there 
are more gluino-degrees of freedom to be produced. 

• Gluino-Squark pair production (gq) 

The gluino-squark pair production process is al- 
ways initiated from a qg or qg. This process is 
dominated by the t-channel exchange of a squark 
or a gluino, but also the s-channel quark mediated 
diagram is available. The cross section of the qg 
process in the Dirac model is identical to the Ma- 
jorana case. A full list of the relevant Feynman 
diagrams and the related cross section formulae for 
all the processes described here can be found for 
example in [5]. 

Let us then look at the total cross section of pair pro- 
duction of supersymmetric particles as a function of the 
squark and gluino masses. Obviously, the cross section of 
these pair production processes is different in the M = 1 
and N = 2 cases, and the total production cross sec- 
tion of supersymmetric particles behaves differently as 
a function of the squark and gluino masses. In the low 
gluino mass region, where gluino pair production dom- 
inates, the total cross section is enhanced in the Dirac 
case with respect to the Majorana case. On the other 
hand, in the low squark mass region the effect is oppo- 
site. This is illustrated in Fig. [TJ where the ratio of the 
production cross sections of supersymmetric particles in 
the MSSM and in the Dirac model is plotted as the func- 
tion of the gluino mass. The squark mass is set to ^, 1 
or 2 times the gluino mass. As explained above, for a 
small squark mass, niq = \rrig the ratio grows rapidly 
for increasing gluino mass. This means that the over- 
all production cross section of supersymmetric particles 
is heavily suppressed in that region of parameter space 



in the Dirac model. This will become important when 
looking at the current bounds from the LHC, Sec. \V\ 



Trig = rrig/2 




500 1000 1500 2000 

m s (GeV) 

FIG. 1: The ratio of cross sections for J\f = 1 and N = 2, 
where the total cross section is given by a = o(qq) + o(qg) + 



IV. SIMULATION 



In the next section we are going to describe how we ob- 
tain bounds on N = 2 SUSY based on the latest ATLAS 



search on multijets and missing energy [TU]. See Sec. V A 
for details. Therefore, we perform a simulation based on 
this experimental analysis. Following their analysis, we 
define a simplified model where the relevant degrees of 
freedom are the first generation squarks u^, Ur, and 
c£r, the lightest neutralino x , assumed to be mostly bino, 
the lightest chargino Xi an d the gluino. The gluino can 
be either a Majorana particle, as in the MSSM, or a Dirac 
particle go, which is equivalently described by the two 
Majorana mass eigenstates g\ and g2 - See Appendix for 
details. We set the masses of the other supersymmetric 
particles to 10 TeV, effectively decoupling them from the 
theory. The bino, which is the lightest supersymmetric 
particle (LSP) and thus stable, is assumed to have a mass 
of 100 GeV. The chargino mass is set to 150 GeV, but 
it is only relevant for Sec. |VI| We then explore different 
configurations for the squark and gluino masses. 

We use Madgraph [11] as the first step to generate 
parton level Monte-Carlo samples. We introduced in 
Madgraph the Dirac gluino model using an interface 
with Feynrules [12]. We also use the interface with 
BRIDGE [13], to handle long decay chains. The samples 
are generated at 7 TeV, with the baseline cuts used by the 
ATLAS analysis [TD]. We then shower the samples with 
PYTHIAv6.4 [14], and use a ROOT-based macro to im- 
plement some of the detector effects from a detector fast 
simulator called ATLFAST [15 ]. Those detector effects 
include jet clustering and energy-momentum smearing. 
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V. THE EXCLUSION LIMIT OF SQUARK AND 
GLUINO MASSES 

A. ATLAS analysis 

In this section we briefly describe the analysis done 
by the ATLAS collaboration in the multijets and miss- 
ing transverse momentum channels [10] with 1.04 fb -1 of 
data. The selection is purely hadronic and a lepton veto 
is applied. The exclusion limits are shown in a (mg, rriq)- 
plane, where all other supersymmetric particles except 
the neutralino are decoupled. 

Their strategy is to separate the signal in five bins of 
different number of jets and missing energy cuts. In all 
bins, the leading jet has px > 130 GeV, and jets with 
Pt > 40 GeV are used for the counting of number of 
jets, rij. The variable M e ff is built by summing missing 
energy and the jets characterizing the different signal re- 
gions. Additional cuts on the separation between missing 
momentum and jets, and on the ratio of missing energy 
to M e ff are applied. The signal regions correspond to 
ri j 2 and 3, two more regions with rij 4 and two 
different cuts on M e ff, and a final region with at least 
four jets with p T ^ 80 GeV and M e// > 1100 GeV. 



B. The exclusion limits for Dirac gluinos 
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FIG. 2: The exclusion limit in the (m,g, m^)-plane in the Af = 
1 (solid line) and in the Af = 2 model (dashed line). 



same gluino mass, and varied the squark mass to get the 
same number of events in the exclusion bin. The points 
in the (nig, m,j)-plane obtained by this method were then 
connected to form the exclusion limit in the Dirac model. 



As explained in section III the overall production cross 
section of supersymmetric particles behaves differently in 
the Af = 1 and Af = 1 cases. Since searches are based 
on the assumption that SUSY is realized as Af = 1 , 
this feature affects the sensitivity of the searches for su- 
persymmetry and the exclusion limits in the (m,g,mq)- 
plane. The exclusion limit in the MSSM (Af = 1 ), as 
obtained by ATLAS [TU], is shown by the solid line in 
figure [2] For the Dirac case this changes, because the en- 
hanced production cross section of gluino pairs increases 
the sensitivity of the search in the low gluino mass region, 
and the suppressed squark pair production cross section 
decreases the sensitivity in the low squark mass region. 
Therefore the exclusion limit becomes tilted, favoring a 
heavy gluino and leaving more room for a light squark. 
The exclusion limit in the Dirac case is shown by the 
dashed line in Fig. [2j 

We obtain this limit by the following procedure: The 
ATLAS search uses five overlapping bins, divided by 
number of jets, the effective mass of the event and the 
transverse momentum of the jets. For each point in 
the (rrig, m^)-plane, the bin with the best sensitivity is 
used to obtain the exclusion limit. For details, see [TU] , 
We produced a sample of events with MadGraph us- 
ing the same model assumptions as the ATLAS search 
(see Sec. V AJ ). We then run a simulation as described in 
Sec. IV applied the same cuts as the ATLAS analysis 
and selected the bin with the highest number of events 
as the bin used for the exclusion. For each point, we then 
produced a sample of events in the Dirac model with the 



VI. DISTINGUISHING BETWEEN DIRAC AND 
MAJORANA GLUINOS 

In the previous section we saw that if the gluino is a 
Dirac-particle, the exclusion limits set by ATLAS and 
CMS do not apply as such. There is a notable reduction 
in the exclusion limit on the squark mass, especially if 
the gluino is relatively heavy, of the order of 2 TeV 1 . 
Thus, we have shown that in Af = 2 SUSY it is still 
possible that first and second generation squarks exist 
at or below the TeV-scale. In this section we turn our 
focus to the question of identifying the Dirac or Majorana 
nature of the gluino, provided experimental evidence of 
the existence of supersymmetric particles is found 2 . 



In section III we described the three parton level fi- 
nal states relevant for the production of supersymmetric 
particles. The squarks and gluinos then cascade decay 
to Standard Model particles and the LSP. If one im- 
poses a lepton veto, the final states we are looking at 
are multiple jets and missing transverse momentum such 
as qq -> 2j + tf T , qg -> 3j + # T and gg -> 4j + $ T . 

However, by imposing a lepton veto we miss other de- 



1 Note that, as the gluino becomes heavier, the difference between 
the Majorana and Dirac cases is reduced. But at moderate val- 
ues of the gluino mass around 2 TeV, the effect of a t-channel 
exchange of the gluino is still dominant. 

2 For a different proposal of distinguishing J\f = 1 and J\f = 2 
SUSY at colliders, see Ref. [9]. 
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cay chains which tell us about the chirality of the squark. 
Indeed, while the right-handed squarks decay almost ex- 
clusively to a quark and the LSP, the left-handed squarks 
have a large branching ratio to — > QXi > where the 
chargino then decays into a W plus the LSP. The W 
can then decay hadronically or leptonically. Hence the 
detector level final state depends on the chirality of the 
produced squarks, with left-handed squarks producing a 
higher jet and/or lepton multiplicity. 

This feature can be used to distinguish between Af = 
1 and TV = 2 . In the MSSM (Af = 1 case) the two 
jets final state mostly originates from the production of 
two right-handed squarks, but in the Dirac scenario this 
process is forbidden. Therefore the two jets final state 
is only accessible via the production of a right handed 
squark and a right handed antisquark, which requires an 
antiquark in the initial state and is thus suppressed by 
the small antiquark PDFs, or through a subdominant 
decay of a left handed squark in the qLq'n final state. 
Therefore the two jets final state is suppressed in the 
Dirac case compared to the MSSM. On the other hand 
the gluino pair production process, which leads to final 
states with a high jet multiplicity, is enhanced in the 
Dirac model. Thus we get more events with rij > 4 
compared to the MSSM. The overall effect is that the 
Dirac case tends to produce more events with a high jet 
multiplicity and less events with a low multiplicity than 
the Majorana case. This is also apparent in table|TJ which 
shows some of the results of our simulation. The number 
shown for each bin is the ratio of production cross section 
of squarks and gluinos times the acceptance of each bin 
for the Af = 1 and Af = 2 cases. The bin with the 
highest acceptance for each value of squark and gluino 
masses is highlighted with a bolded (underlined) font in 
the Af =1 (Af = 2 ) case. The cross sections aree based 
on a tree level calculation performed with MadGraph, 
and the overall scale would have to be adjusted by a 
computation at NLO, but this iG factor is common to 
Af = 1 and Af = 2 and therefore it factors out in the 
ratio of events. One can see a clear tendency towards the 
higher multiplicity bins in the Dirac model. The ratio 
of the number of events in the 2j + TpT final state and 
the number of events with four or more jets plus missing 
transvere momentum can then be used to separate the 
Dirac scenario from MSSM. A low number for this ratio 
is characteristic to the Dirac model. 

To obtain the best current limits, we have followed the 
ATLAS search which assumed a lepton veto. However, 
for the purpose of distinguishing between Af — 1 and 
Af = 2 , the number of events with one or more leptons 
compared to the number of events without leptons is a 
useful observable. If we allow for one or more leptons 
in the two jets plus missing transverse momentum bin, 
the qLQR final state may end up in this bin via the decay 
Ql — > qxf ~ * ql ±l/ ix - Since the qnq' R partonic level fi- 
nal state is not accessible in the Dirac model, this decay 
chain makes up for a notable part of the two jets events. 
Therefore the ratio of events in the two jets bin that con- 
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TABLE I: The ratio of the number of events in A/" = 1 /A/" = 2 
SUSY models, for the bins used in the ATLAS search for 
hadronic events and missing transverse momentum |10j . 



tain one or more leptons to the number of events in that 
bin that do not contain isolated leptons is high in the 
Dirac model. In Af — 1 SUSY the majority of events in 
the two jets bin is produced via the qBq'n partonic level 
final state. Therefore allowing for one or more leptons 
does not add significantly to the total number of events 
in that bin, since the leptonic decay of the left handed 
squark in the qL<j_R final state is a subdominant process 
compared to the hadronic decay. Thus a good separation 
between the models is achieved by comparing the ratio 
of events in the two jets plus missing transverse momen- 
tum bin that do or do not contain one or more leptons. 
This ratio is higher in the Dirac case compared to the 
Majorana case. 

The best separation between Af = 1 and Af = 2 is 
achieved by combining these two observables. This is 
schematically shown in figure [3] where the ratio of events 
in the two jets and four jets bins is plotted in the x-sods 
and the ratio of events with or without leptons in the two 
jets bin is shown in the y-axis. We generated event sam- 
ples with different mass configurations, but all obeying 
rriq < trig (see discussion below for the vtiq > nig case). 
We then placed the resulting point of each mass config- 
uration on the (x, y)-plane. The points corresponding to 
Af = 1 end up on the lower right corner of the plane as 
expected, and the points in Af — 2 populate the upper 
left corner of the plane. When the gluino mass becomes 
very heavy, however, the gluino starts to decouple from 
the theory. Then the squark pair production process be- 
gins to be dominated by the s-channel gluon exchange 
diagram, which contains no information of the Dirac- or 
Major ana-nature of the gluino. Therefore the points cor- 
responding to the Dirac model begin to approach the 
MSSM part of the (x, y)-plane in the limit where the 
gluino mass is very heavy. This is also illustrated in table 
PH where the values of these observables is shown for two 
points. One of the points is in the area where the gluino 
is not too heavy and the separation between Af = 2 and 
Af = 1 is good, and the other one shows how the separa- 
tion gets worse as the gluino mass grows and the squark 
mass gets small. We conclude that, except for the case 
where the gluino is extremely heavy and thus practically 
decoupled from the theory, the two observables plotted 
in the x and y axes of the plane can be effectively used to 
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TABLE II: The relevant observables used for the identification 
of the Dirac nature of the gluino: the number of events with 
two jets versus the number of event with four or more jets, 
and the number of events with one or more leptons and two 
jets versus two jets and no leptons, in the Af = 1 and Af = 2 
cases for two points in the (mg,mq) plane. 



identify if the gluino is a Dirac- or a Majorana-particle. 

The above picture is less clear if the gluino is lighter 
than the light squarks, because then all the squarks de- 
cay dominantly to a quark and a gluino, which leads to 
higher jet multiplicities also in the Af = 1 case. The lep- 
ton vs no-lepton discriminant is still a good discriminant 
between the two scenarios, but the ratio of two jets finals 
states and four or more jets final states is very small in 
both models in this region of the parameter space, and 
one would have to look into observables with two more 
jets, namely the ratio of four to six jets versus the ratio 
of lepton and no lepton in the four jet final state. 

Let us finish by mentioning the option of pseudo-Dirac 
gluinos, i.e., Dirac gluinos with a small Majorana admix- 
ture. A small splitting between the two Majorana states 
in the Dirac gluino would not affect the results in this 
section, as they are based on multijet final states with 
high-pT requirements. 
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FIG. 3: A schematic illustration of the two observables used 
for distinguishing between the Dirac model and the MSSM. 
On the a;-axis is the ratio of the number of events with two 
jets and the number of events with four or more jets. On the 
i/-axis is the ratio of events with or without leptons in the two 
jets bin. The two areas shown are obtained by varying the 
squark and gluino masses in the two models. When the gluino 
is much heavier than the squark, both Af = 1 and Af = 2 cases 
become indistinguishable. 



tween Af = 1 and Af = 2 supersymmetry. Since nat- 
ural models of supersymmetry are realizations of either 
Af = 1 or Af = 2 in the gauge sector, this discriminant 
is effectively counting the number of supersymmetries in 
Nature. 



VII. CONCLUSIONS 

In this paper we have shown that the current bounds 
on Supersymmetry are drastically altered when one con- 
siders the possibility that Supersymmetry could be re- 
alized as an Af = 2 version in the gauge sector. This 
is a key question nowadays, as the LHC is pushing the 
bounds on Supersymmetry to uncomfortable areas for 
natural solutions of the hierarchy problem. 

We have taken the best current limits from ATLAS 
and redone their analysis -within the tools available to 
theorists- to estimate the limit on Af = 2 Supersymme- 
try. We find that, in Af = 2 supersymmetry, squarks 
as light as few hundreds of GeVs are still allowed. This 
study illustrates the following important point, that ex- 
perimental studies of Supersymmetry are remarkably de- 
pendent on the assumption about the number of super- 
symmetries. 

We have gone a step ahead, and found a suitable com- 
bination of observables which, provided an excess in mul- 
tijet events is found, can be used as a discriminant be- 
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Appendix 

Here we will show in detail how the Af = 2 gluino sec- 
tor can be expressed in terms of one gluino which has a 
Dirac mass term. We will start by defining some nota- 
tion and conventions, then we show how the Lagrangian 
is transformend from the flavor basis to the mass basis, 
expressed in terms of the two Majorana mass eigenstates. 
Finally we show that this is equivalent to the Lagrangian 
of one Dirac particle. 
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We use the chiral basis, and the notation in [IB]. We 
define two-component Weyl spinors with dotted and un- 
dotted indices: r\ a and r/ a , as complex, anticommuting 
objects with a = 1, 2 and a — 1, 2. The spinors with 
dotted and undotted indices are related via Hermitian 
conjugation: 

Vl ■= (riaV = tfh, ^ tQ = (^) t - (1) 

The spinor indices are raised and lowered using the 
antisymmetric symbol e 12 = — e 21 = — £12 = £21 = 1, 
eu= e» = 0: 



(2) 



and similarly for the spinors with dotted indices. In the 
following, we will suppress contracted indices. 

A four-component spinor is written in terms of the 
two-component spinors as 



* = 



Va 



(3) 



and thus its conjugate spinor is 

* = *V = (vl X a ) 7° = (x a vl) ■ (4) 

In our conventions, the left and right handed fields are 
then simply the upper and lower components of the four- 
component field: 



• 'i 1 • * R =p R * = [ JJ. (5) 







For SM Dirac fields, such as a quark field Q, we therefore 
write 



and 



The charge conjugate of a 4-component spinor is 

— (5$ -(*)■« 

Majorana particles must obey ^ c = 'i', implying the Ma- 
jorana condition X = V- 

Let us move into the gluino mixing. Assuming we have 
two Majorana gluinos: 



h ;/4 > 



(9) 



The mass term in the general case is 

An = ( m 39'g' + m 399 + m D{g'g + 99')) ■ (10) 

We will be interested in the Dirac limit, i.e. m 3 = 7713 = 
0. In terms of the two-component spinors the mass La- 
grangian is then 

C m = --m D (g'g + gg') = -m D {r{r] + r?'V)- (H) 
The mass eigenstates are given by 



9i 



^2^ + ^y ^ 



^2 Ut« _ 



(12) 



(13) 



These are Majorana particles, as g\ — g\ and g\ = gi- 
One can easily check that the mass Lagrangian is indeed 
diagonal in this basis. 

We will now examine the original Lagrangian term by 
term, transforming into the mass eigenbasis. First, the 
kinetic term: 



c K = ^m + g'09') 



= \{X^9^ + X f ^ X + X'^d^ + X ,f ^ X ') 



= ^{h09i + 92^92)- 

Similarly, the gauge term is 

/We = -Wsf ljk Cg l J ~g k +g H 3 g' k ) 
= -i9sf ij \m j yi + &@ i ~9l)- 



(14) 



(15) 



The degrees of freedom present in gi and gi can be 
combined into one Dirac spinor gu, given by 



1 

71 v 



9d = -^(51 -ifa) = ( ^ 



(16) 



In terms of the Dirac gluino, the mass term is 

Cm = -m D (r)v' + vW ) = -rn D g D g D , (17) 
as expected. The gauge term and the Kinetic term are 



C 



gauge 



(18) 



Ck = 2 (gD0gD + ~g c D $g c D ) 



The trilinear quark-squark-gluino coupling term is 

C% = ~9s {Ql{)d Ql - Qn.g c D qR) + h.c. . (19) 
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Thus we have seen that the gluino sector is equivalently 
described either by two degenerate Majorana particles 
c/i and §2 or by one Dirac particle <?£>. To validate our 
Monte Carlo implementation, we implemented the M = 2 
model twice, using both descriptions, and checked that 



the numerical results agree between these models. We 
also checked that decoupling one of the Majorana gluinos 
reduces the TV = 2 model back to MSSM, and our nu- 
merical results confirm that. 
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